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We study oscillatory universes within the context of Loop Quantum Cosmology. We make a com- 
parative study of flat and positively curved universes sourced by scalar fields with either positive 
or negative potentials. We investigate how oscillating universes can set the initial conditions for 
successful slow-roll inflation, while ensuring that the semi-classical bounds are satisfied. We observe 
rich oscillatory dynamics with negative potentials, although it is difficult to respect the semi-classical 
' bounds in models of this type. 
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The study of cyclic/oscillatory universes has a long history. One of the primary motivations for such a study was to 
circumvent the need for initial conditions in cosmology. This was, however, shown to be extremely difficult to achieve 
_ | within the context of general relativity (GR), without encountering singularities. 

Recent developments in String/M-theory inspired braneworld models have revived hopes for the possibility of 
constructing such universes (see e.g. Ref. One such attempt is the so called cyclic/ekpyrotic scenario (see Ref. 

y—{ , and references therein) . Despite the attractiveness of this possibility, a successful (non-perturbative) treatment of 
■ bounces within the context of M-Theory is still lacking. 

More recently such universes have been studied in the context of loop quantum gravity, which at present is the 
CI main background independent and non-perturbative candidate for a quantum theory of gravity (see for example Refs. 

and 0). To date, cosmological applications have focused on minisuperspace/midisuperspace settings with a finite 
I . number of degrees of freedom. This program is referred to as Loop Quantum Cosmology (LQC). In this context the 
evolution of the universe can be divided into 3 distinct phases: an initial high energy and high curvature quantum 
phase, described by a (or a set of) difference equation(s); an intermediate semi-classical phase, with continuous 
evolution equations modified due to non-perturbative quantization effects; and finally, a classical phase, where the 
usual continuous cosmological equations are recovered. It has recently been shown that such semi-classical effects can 
result in bouncing FRW universes which avoid singularities |f| . This has in turn led to the possibility of non-singular 
cyclic universes in LQC settings 0, Q • 

In this paper we make a comparative study of recent results concerning both spatially flat ||| , and positively curved 
oscillating universes [(|, with both positive and negative potentials. We extend these results in a number of areas, 
in particular, in the case of models with negative potentials. This latter class of models is of direct relevance to the 
5^ 1 ekpyrotic/cyclic models considered recently. 
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II. EFFECTIVE FIELD EQUATIONS IN LOOP QUANTUM GRAVITY 

The semi-classical phase in isotropic LQC refers to the regime where the scale factor lies in the range < a < a*, 
where = y/~flp\, a* = \J 7j/3Zpi, 7 = ln2/-\/37r « 0.13 and j is a quantization parameter which must take half 
integer values. Below the scale a^, the discrete nature of space-time becomes important, whereas the standard classical 
cosmology is recovered above a*. The parameter j therefore sets the effective quantum gravity scale. The modified 
Fricdmann equation is given by 



H = \ - \ = 
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where k takes the values or 1 for a flat or positively curved universe, respectively, and the quantum correction factor 
D(q) is defined by 



D(q) = (±\ q^ {7 [{q + l) 11 / 4 - \q - l| n / 4 ] - llq [(q + l) 7 / 4 - Sgn(q l)\q 1| 7 / 4 ] } 
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with q = (a/a*) 2 . As the universe evolves through the semi-classical phase, this function varies as D oc a 15 for a <C a*, 
has a global maximum at a«a„ and falls monotonically to D = 1 for a > a* . 
The scalar field equation has the form 

where a prime denotes differentiation with respect to the scalar field. Differentiating Eq. and substituting from 
Eq. J2I then gives 

We can immediately see two interesting features. First, the correction to the scalar field equation © causes the 
frictional/anti-frictional 4> term to change sign when dlnD/dlna passes through 3. Secondly, an expanding universe 
naturally undergoes a period of super-inflationary expansion (H > 0) pi when d In D/d In a > 6 (assuming the 
curvature term to be either zero or negligible). 

An interesting consequence of these features is that both a collapsing positively curved universe 0, , and a flat 
universe sourced by a scalar field with a negative potential can undergo a non-singular bounce. We can most 

clearly see how the bounce arises by reformulating Eqs. and J2J) in the standard form of the Einstein field equations 
sourced by a perfect fluid [||: 

Tj2 _ 87T& k 
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PcS = -3H (poff + Pes) , (6) 

where 

Pe s^- 2 ^ + V, (7) 
10 2 / ldlnZA 

55 slri 1 -^)-^ (8) 

define the effective energy density and pressure of the fluid, respectively. The effective equation of state, w = p c g/ p e e, 
is given by 

20 2 / IdlnlA 
u; = -l + - 1 . 9 

<j) 2 + 2DV V 6 dhxa J y ' 

Hence, the LQC corrections to the cosmic dynamics can be entirely parametrized in terms of the equation of 
state. When dlnD/dlna > 6, the fluid represents 'phantom' matter (w < —1) that violates the null energy condition 
{peS+PeS > 0), independent of the form of the potential. During a collapse the energy density starts to decrease when 
this condition is met, and in a positively curved universe, the density term in the Friedmann equation is eventually 
balanced by the growing curvature term, thereby leading to a bounce. In the case of a flat universe sourced by a scalar 
field with a negative potential, the decreasing energy density instantaneously vanishes as the positive kinetic energy 
component is balanced by the negative potential. In both cases, H tending to zero necessarily leads to a bounce 
because H > 0. 

In the following sections we shall study the question of initial conditions for inflation in both flat and positively 
curved universes, with both positive and negative potentials and we will investigate whether other interesting scenarios 
involving negative potentials can be developed. 



III. SCALAR FIELDS WITH POSITIVE POTENTIALS 



In this section we consider a semi-classical universe sourced by a scalar field self-interacting through a positive 
potential. We are particularly interested in how the dynamics of the semi-classical region can set the initial conditions 
for slow-roll inflation if the field is situated at or near the minimum of its potential. Our discussion applies for an 
arbitrary potential apart from the weak assumptions that it has a global minimum at V^n n (0) = and is a positive- 
definite and monotonically varying function when 4> ^ such that V" > 0. We shall also assume that the field is 
initially located at the minimum of its potential. For definiteness, we focus in this section on a massive field with a 
quadratic potential, V((f>) = m 2 2 /2. 
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Flat case 



For a flat universe the anti-friction mechanism simply accelerates the scalar field up the potential [ME!- A major 
problem with this mechanism is that during the anti-friction era the energy density of the field increases. On the 
other hand, the Hubble length must remain larger than the limiting value of the scale factor consistent with the 
semi-classical regime, i.e., \H\ai < 1. (This is roughly equivalent to the requirement that energy scales remain below 
the Planck scale and is referred to as the Hubble bound). Since this length decreases during the period of anti- friction, 
we need to determine how far the field can be moved before this condition is violated. To achieve 60 e-folds of inflation 
for a quadratic potential, V = m 2 (f) 2 /2, the field must be displaced by at least 3lp^ from the minimum of its potential 

at the onset of inflation. It was recently found that given an initial value <^i n it, the field could only be moved as 
far as 2.4/pj 1 or less (for sufficiently large j) without violating the Hubble bound. 

The behavior of the field can be understood by approximating the evolution into two epochs. In the first, between 
ainit and the end of the semi-classical regime as, the asymptotic form, -Dapprox = (12/7) 6 (a/a*) 15 , is used. In the 



second with a > ag, we assume D — 1. Here as is defined as the value of the scale factor where D a 



first reaches 



unity. The field is assumed to be massless up to the turning point in its evolution, which is estimated to occur when 
m 2 (j> 2 /2 w <pi, where the subscript t stands for "turning point". Hence, neglecting the potential in the Friedmann 



equation QJ, we have (dcp/dlna) 2 
obtain 



6-D/87rZp[. Integrating for <p, and using the approximation above for D, we then 
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The explicit dependence on at can be removed by using the expression 



as 
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which itself arises from using the first integral 



4> = 0init (^^) 
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in both epochs, along with the form for -Dapprox- Finally, applying the condition for the turning point, we arrive at 
the estimate 



h exp 



V2 / 7 
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We note that for a given initial value q- m i t — (ai n it/ a *) 2 i the value of 4> when inflation occurs is independent of j 
and proportional to </>; n it. This approximation was compared to numerical results @, and found to be in very good 
agreement. We emphasize that this approximation only works for initial conditions such that ainit *C a* because the 
approximation for D involves that assumption. 



B. Positively curved case - massless scalar field 



Before considering the positively curved case with a positive potential, it is illuminating to consider the dynamics 
of such a universe sourced by a massless scalar field. This was discussed in detail in Ref. [(j. We can, however, easily 
explain the form that the dynamics will take by the following simple argument. We know from the above discussion 
that a collapsing positively curved universe will always bounce into a new expanding phase no matter what form 
the potential takes, and consequently, this bounce will still occur even for V — 0. We also know that an expanding, 
positively-curved, classical cosmology sourced by a massless scalar field will rapidly undergo a re-collapse since the 
energy density scales as p oc a~ 6 , and therefore redshifts much more quickly than the curvature term in Eq. (J3J). After 
the bounce, the universe expands and once it has entered the classical phase a re-collapse will inevitably occur. The 
process then repeats itself indefinitely. 

Since the effective equation of state (0 is independent of the field's kinetic energy when V = 0, we see that identical 
cycles must occur. Furthermore, Eq. (|12fl implies that the field's kinetic energy never vanishes during a given cycle 
since the scale factor always remains finite. The value of the field therefore increases or decreases monotonically with 
time. 
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C. Positively curved case — massive scalar field 



To understand how the cyclic dynamics is modified by the presence of a potential, we recall that the equation of 
state can change from cycle to cycle depending on the field's position on the potential. Assuming that the potential 
is initially unimportant, the dynamics is qualitatively similar to that described in section 3.2. This implies that the 
kinetic energy of the field never vanishes and a re-collapse always occurs. However, since the field moves monotonically 
away from the minimum, the potential becomes progressively more important with each completed cycle and it must 
therefore become dynamically significant after a finite number of cycles have been completed. In general, a re-collapse 
occurs if the matter redshifts more rapidly than the curvature term. This requires the strong energy condition to 
hold, w > —1/3, which in turn implies V < <fi 2 . As the field moves up the potential, w is pushed towards —1 and it 
becomes progressively harder to satisfy the strong energy condition. Eventually, therefore, a cycle is reached where 
V > (p 2 and the condition w > —1/3 is violated. This leads to an epoch of slow-roll inflationary expansion, with the 
field slowing down, reaching a point of maximum displacement and rolling back down the potential. For some initial 
conditions and choices of parameters, this may occur during the first cycle, in which case the scenario discussed for 
the flat model is recovered. 

Thus far, we have discussed the cyclic dynamics for a given value of the quantization parameter, j. In spatially 
flat models, we saw that the maximum value attained by the field as it moves up the potential is independent of 
j given an initial ratio ai n it/a» and 0i n it- For positively-curved models, on the other hand, the field moves further 
up the potential before the onset of slow- roll inflation for smaller values of j. This behavior can be quantitatively 
understood by assuming that the total energy at the point of collapse is conserved through the cycles and that the 
potential is dynamically insignificant during the first cycle. The massless case, therefore, gives us a good estimate of 
the total energy, with the evolution of this energy being given by Eq. I|12p. Now, assuming the universe to be classical 
at <frt , using the condition V t ~ <p 2 to estimate the point at which inflation starts, and the relation 87W P1 pi/3 k, 1/a 2 
to eliminate a t , one obtains 



2 J - / 87rZpA Anit 



3/2 
9init 



(14) 



This demonstrates that for the quadratic potential, with a given initial a; n it/a*, the value of V when inflation occurs 
is inversely proportional to a 3 and hence inversely proportional to j 3 / 2 . It also shows an inverse dependence on <^>i n it - 
It is worth pointing out that in deriving Eq. 114f) . unlike in the flat case, we have not used any approximation to 
the general form of the correction factor D(q), and consequently this expression is still a good estimate of 4>t when 

^init ^* • 

As in the flat model, it is important to ensure that the Hubble bound is not violated. However, due to the inverse 
dependence of (f>t on (/>; n it, this proves to be a much weaker constraint in the positively curved models. The field is 
moved up the potential in a series of small 'kicks' of anti-friction, and can easily be moved sufficiently to drive 60 
e- folds of slow- roll inflation in the case of a quadratic potential |H| • 

In Fig.Q w e illustrate the dependence of the maximum value attained by the field before turnaround on the value of 
the quantization parameter j and the ratio a; n it/«* for a given initial value of the kinetic energy </>i n it. The horizontal 
lines represent the estimated value of <pt extracted from Eq. <]13[l and the lines with negative slope represent the 
estimate Ijl4(l . Circles and triangles correspond to the actual values obtained by numerically integrating the equations 
of motion. When many oscillations occur before the onset of inflation, Eq. (|14|) yields a good estimate of 4>t- However, 
if the field is moved far enough in the first period of anti-friction for the strong energy condition to be violated, no 
oscillations occur and the situation is equivalent to the flat model. In this case, Eq. I|13fl gives a good estimate of </>t. 
The points where the lines cross in Fig. ^ show the transitions between the two regimes. 

The numerical results in Fig. ^show a stepping behavior, with a range of j leading to roughly the same value of 
</>t, before there is a drop to the next value. This becomes more pronounced at higher j. This occurs because our 
approximation relies on two estimates: firstly, that V ~ 4> 2 at the onset of inflation, and secondly, that 87rZ|,jp f /3 « 
1/a 2 when the field reaches its maximal value. Clearly, the two approximations can not be satisfied simultaneously, 
since the first implies a violation of the strong energy condition, whereas the second represents a turning point in the 
expansion. In effect, since we assume that the potential becomes dynamically significant only during the last cycle 
0, we are carrying over the second condition from the previous cycle. Since 87rZpjp t /3 > 1/a 2 for the cycle in which 
inflation occurs, the second condition underestimates the amount of energy present and consequently underestimates 
the change in the value of the field before it turns around. A range of values of j will give rise to the same number 
of cycles: for the lowest value of j in a given range the condition 87rZ P1 pt/3 « 1/a 2 is a good approximation, since a 
re-collapse is only just avoided, but for subsequent values the underestimate becomes progressively more severe. This 
gives rise to the steps which become more evident towards larger j . This can be understood by recalling from Eq. i|14|) 
that 4>t oc j -3 / 4 , which implies Aj cx j 7 / 4 A</>, thus showing that the steps in j are wider for larger j. 
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FIG. 1: Illustrating the dependence of the maximum value of the scalar field, 4>t, on the quantization parameter j for the 
quadratic potential V(<j>) = m 2 (j> 2 /2 with m = 10 — 6 Zp 1 1 . The horizontal lines represent the estimated value of cf>t derived from 
Eq. 113|l and the tilted lines represent the estimate l|14|l . The solid lines correspond to the initial conditions ai a it/a* = 0.7 and 
the dashed lines to a[ n it/a t = 0.5. Circles and triangles represent the actual values obtained by numerically integrating the 
equations of motion for ai n it/a* = 0.7 and a- m n/a» = 0.5. 

IV. SCALAR FIELDS WITH NEGATIVE POTENTIALS 

In the spatially flat setting, we have seen that a bounce can occur when the universe is sourced by a scalar field 
self-interacting through a negative potential. It is important to investigate whether interesting cosmological scenarios 
can be developed using negative potentials. For example, a negative exponential potential plays a key role in the 
cyclic/ekpyrotic scenario 2]. Alternatively, bouncing cosmologies in which the maximal size of the universe increases 
with each successive bounce could offer a possible resolution of the flatness and horizon problems of standard, big 
bang cosmology (see Ref. ^3 an d references therein). 

Here we look at the general dynamics of a cosmic bounce with a negative potential. During the collapsing phase, 
once the scale factor falls sufficiently for the condition d\nD/d\na > 3 to be met, the kinetic energy of the scalar 
field decreases due to the friction in Eq. @. This is eventually balanced by the negative potential energy and leads 
to a bounce. Since the balance is achieved before the kinetic energy reaches zero, the field continues to move in 
the same direction during the bounce. This is similar to the dynamics of a positively curved model with a positive 
potential. After the bounce, a super-inflationary era ensues, and the standard scenario is recovered once the scale 
factor exceeds a*. The field subsequently decelerates, and in the case where the potential is negative-definite, a 
re-collapse is initiated when the kinetic energy is again balanced by the potential. The qualitative behavior discussed 
above is then repeated indefinitely with the field always moving in the same direction. The behavior, however, is 
not symmetric since it responds to the region of the potential where the field is positioned. Consider, for example, 
the case of an exponential potential, V cx — exp(— a<f>): if the field is moving up the potential (towards more positive 
values), then the maximum size attained by the universe increases with each cycle whereas the minimal size decreases. 
This occurs because the universe has progressively more time to evolve before the kinetic energy is canceled by the 
decreasing negative potential term, and as the potential asymptotes to zero, the maximal/minimal size of the universe 
stabilizes. Conversely, if the field is moving towards more negative values, the reverse is true. 

One question that naturally arises is whether this mechanism for realizing a non-singular bounce may be employed 
within the context of the cyclic/ekpyrotic scenario? This question has been asked previously [7(, and the discussion 
given here is in agreement with the previous study. There are three major issues that need to be addressed. Firstly, 
in the cyclic/ekpyrotic scenario, the field is supposed to change direction at the bounce, but as we have discussed, this 
is difficult to achieve at the semi-classical level with the mechanism under discussion, although this problem might be 
alleviated by introducing a steep positive section to the potential. (This would cause the field to slow down and then 
reverse its direction of motion). 

The second problem is more severe and is directly related to Hubble bound Hdi < 1 discussed above. During 
a classical collapse, the kinetic energy of the field increases due to the anti-friction in the field equation J3J). The 
kinetic energy can start to fall only when the condition d In D/d In a > 3 is attained. Consistency at the semi-classical 
level requires that the Hubble bound must not be violated before this point is reached. (This problem of Planck scale 
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energies has been raised previously [Uj-) In the case of a negative exponential potential, there exists a scaling solution 
in which the kinetic and potential energies of the field vary in proportion to one another. In principle, therefore, the 
total energy density may remain below the Planck scale. For the cyclic/ekpyrotic scenario, however, the potential has 
an exponential form only over a finite region or parameter space, and once the field evolves through this region, the 
cosmic dynamics moves away from the scaling solution. This implies that the kinetic energy is no longer balanced by 
the negative potential, and the total energy rapidly exceeds that of the Planck scale. 

For a collapsing universe sourced by a massless scalar field, it is possible to estimate when the bounce must occur in 
order for the Hubble bound to remain satisfied throughout the evolution. In particular, we can derive a lower limit on 
the value of the ambiguity parameter, j. If the collapse starts well into the classical regime, it follows that 7Ji n it = 0, 
D = 1, and that the initial kinetic energy of the field is 0? nit = & / %-Klp y a 2 n{t . The subsequent evolution of the field is 
determined from the relation <j) 2 = </>f n ; t (ainit/a) 6 and, after substituting this expression into the Friedmann equation 
(|T|) and imposing the consistency condition Ha,i < 1, we require that 

J>-(^-)\ (15) 

where a* has been employed as the point at which the energy density becomes maximized. Clearly, if the size of the 
universe is large compared with a* at the start of the collapse, j must be very large indeed if a bounce is to occur 
within the validity of the semi-classical regime. Although a realistic collapse will involve a massive (or more generally 
a self-interacting) scalar field, the above estimate highlights the difficulties that arise: even if the field equations 
generically lead to a bouncing cosmology, the bounce may occur in a regime where the equations are no longer valid. 

Finally, a third question that arises in developing a consistent bounce is also related to the size of the universe. For 
a consistent, semi-classical treatment, it is important that the bounce does not occur after the universe has collapsed 
beyond the minimal scale consistent with the semi-classical region, a*, since the semi-classical Friedmann equations 
are no longer valid below this scale. In particular, this is relevant to any oscillating scenario involving a negative 
exponential potential in a flat universe where the maximal size of the universe increases with each successive cycle: as 
the value of the field increases with each cycle, the value of the scale factor at the bounce decreases and it is possible 
that it may fall below the critical value <Zj. 

Before concluding, it is also worth considering the case where the potential has a global minimum V m i n , such that 
Knin < 0, with only a region around the minimum that is negative. One example of this class is a potential of the 
form V(4>) = C + m 2 (f> 2 /2 for some negative constant C . In positively-curved universes, the mechanism described in 
section 3.2 enables the field to work its way out of the negative region and continue up the potential until the cycles 
are broken. In Fig. [21 we illustrate the evolution of the scalar field, the scale factor and Hubble ratio for a quadratic 
potential V = C + m 2 (f) 2 /2 in a positively-curved universe. The figure shows the stepwise evolution of the field up the 
potential followed by a short period of inflation between ln(mt) = [1, 1.5]. As the field rolls down the potential, it once 
again reaches the negative section of the potential and this eventually causes the universe to re-collapse. During the 
collapse the Hubble parameter grows and is peaked around a « a*. Below a*, the potential prevents the scale factor 
from becoming singular and the field resumes the same stepwise behavior, but now in the opposite direction. This 
results in a further period of inflation. For spatially flat universes, the picture is similar although the oscillations come 
to an end as soon as the field reaches a positive region of the potential. Unfortunately, in the flat case the resulting 
amount of inflation is typically less that 60-efolds, even for the quadratic potential. (More inflation is possible but 
at the expense of violating the Hubble bound). Furthermore, for such a potential in a flat universe, there is the 
possibility that the potential becomes positive during the collapsing phase, in which case the universe can undergo 
no further bounces and ultimately collapses into a big crunch. 



V. CONCLUDING REMARKS 



We have summarized and extended recent results concerning the effects of loop quantum gravity corrections to 
the cosmological evolution equations. In particular, we have made a comparative study of the ability of oscillatory 
universes to set the initial conditions for successful slow- roll inflation [f| We have considered both flat and 
positively curved models sourced by a scalar field with either a positive or negative self-interaction potential. Negative 
potentials are of interest as they arise naturally in string/M-theory motivated models such as the recently proposed 
ekpyrotic/cyclic scenario Q- We find that the requirement for self-consistency in the semi-classical regime leads to 
strong constraints that any realistic model must satisfy. 
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1.5 

In (mt) 



FIG. 2: Time evolution of <f>, ln(a) and H for a positively-curved universe sourced by a quadratic potential V(4>) = C + m 2 (f> 2 /2, 
Z" 1 , j = 10, a init /a, = 0.5, J7 in it = and C = -lO -11 /" 4 



where m = 10 6 Z P1 1 , j = 10, ainit/i* = 0.5, Hitdt = and C = —10 n £ P 4 . In the figure, the axes are labeled in Planck units. 
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